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Abstract 

In this paper, the singularity formation of classical solutions for the 
compressible Euler equations with general pressure law is considered. 

The gradient blow-up of classical solutions is shown without any small¬ 
ness assumption by the delicate analysis on the decoupled Riccati type 
equations. The proof also relies on a new estimate for the upper bound 
of density. 
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1 Introduction 

We consider the one dimensional compressible Euler equations in the La- 
grangian coordinates: 


/ 

Tt-Ux = 0 , 

< tii -|- Pa; = 0, (1-1) 

where x is the space variable, t is the time variable, u is the velocity, p is the 
density, r = p~^ is the specific volume, p is the pressure, e is the internal 
energy. Due to the second law of thermodynamics, r, p and e are not inde¬ 
pendent, the relation within which is determined by the state equation(c./. 
m)- Normally, another physical quantity entropy S is considered, which 
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formulates the state equation as p = p{t,S). For solutions, the third 
equation of ( 1 . 1 ) is equivalent to the conservation of entropy {c.f. my- 


St = 0. 


( 1 . 2 ) 


Apparently, ( 1 . 2 ) shows that S is just the function of x. 
pressure law we consider in this paper is 


And, the general 


P = p{t,S) =p{t,S{x)). 


(1.3) 


Then the system (1.1) becomes 


Tt-Ux = 0 , 

ut +p{t,S{x))x = 0 . 


(1.4) 


We consider the calssical solution of initial value problem for (1.4) with 
initial data 

T{x,t = 0) = To{x), u{x,t = 0) = uo{x). 


Compressible Euler equations is one of the most important physical mod¬ 
els for systems of hyperbolic conservation laws. It is well known that shock 
waves are typically formed in finite time and the analysis on the system 
is difficult because of the lack of regularity. The singularity formation for 
both the small initial data problem and the large initial data problem has 
long been a very serious issue for the systems of conservation laws. The 
well-posedness theory for systems of hyperbolic conservation laws could be 
found in [DEKinillH]. 

When initial data is small, the singularity formation has been well stud¬ 
ied for decades. Lax [12] proved that singularity forms in finite time for the 
general systems of strictly hyperbolic conservation laws with two unknowns 
with some initial compression. For general systems of conservation laws, 
nn ini (a US] provide fairly complete results for small data. Specihcally, 
these results prove that the shock formation happens in hnite time in any 
truly nonlinear characteristic field if the initial data includes compression. 

However, the large data singularity formation theory has been finally 
established in very recent papers for isentropic Euler equations with 

7 -law pressure {p = and the full compressible Euler equations of 

polytropic ideal gas {p = K^e,^ ! t~"<) ^ where Ki, K2 are positive constants 
and 7 > 1 is the adiabatic gas constant. The key point in proving the finite 
time shock formation for large solution is to have sharp global upper and 
lower bounds of the density. More precisely, if we restrict our consideration 
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on singularity formation for full compressible Euler equations, the uniform 
upper bound of density is needed for any 7 > 1 , but the time-dependent 
lower bound of density is needed only for the most physical case 1 < 7 < 3 
(c./. m)- The uniform upper bound on density for 7 -law pressure case has 
been found by [7] which directs to a resolution of the shock formation when 
7 > 3. The singularity formation problem when 1 < 7 < 3 was finally 
resolved by [ 1 ], in which the authors proved a crucial time-dependent lower 
bound estimate on density lower bound. Later on, the time-dependent lower 
bound of density is improved to its optimal order 0{l/t) in [3]. 

Nevertheless, for the full compressible Euler equations with general pres¬ 
sure law, the singularity formation results for non-isentropic case are still 
not satisfied when the smallness assumption on the initial data is removed. 
In fact, a complete finite time gradient blow-up result has been showed in [1] 
when entropy S' is a given constant. Furthermore, [ 6 ] provides a singularity 
formation result for the non-isentropic general pressure law case. Unfor¬ 
tunately, in [ 6 ], there are still several a priori conditions on the pressure 
function which are not automatically satisfied for the gas dynamics. The 
target of this paper is to establish a better singularity formation result on 
non-isentropic Euler equations without such kind of a priori assumptions. 
The key idea is to establish a uniform upper bound estimate on density, 
which was lack for general pressure law case previously. In this case, the 
lower bound of density is redundant. Our proof relies on the careful study on 
the decoupled Riccati type ordinary differential equations on gradient vari¬ 
ables which was provided in [^. Using our new estimates, we can get the 
constant lower bound on coefficients of the Riccati type equations, and the 
quadratic nonlinearity implies the derivatives must blow-up in finite time. 

Through out this paper, we need to propose the following assumptions on 
the pressure: there exists a positive function m = m{S), positive constants 
A , k > 1, ki, k 2 and k (f = 1, 2, • • • , 8 ) such that, for r G (0, -|-oo). 


(HI) 

(H2) 

(H3) 

(H4) 


Pr < 0, Pt-t > 0, limp(T) = -1-00, lim p{t) = 0, (1.5) 

T^O T^ + OO 


r*l /•+00 

yJ-Prdr = -1-00, J y/-prdT < -|-CO, 


( 1 . 6 ) 


hpc^ < PPtt < hc2 , (1.7) 

2{k - 1){-Pr f > kpprr, (5 A){prT)‘^ - AprPrrr > 0, (1.8) 

m!{x) m'{x) m! , , m! 


k2m 


-p<Pfi< i —Pr < Prfi < — Pt, (1-9) 

kim l^m Ism 


m' m' m' m' 

- p < P^l^l < -p, - Ptt < PrT^l < 7- Ptt- (1.10) 

L^m l^m Ism Ljm 
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Here, the sound speed is 


c = V-Pr{r, S) 


( 1 . 11 ) 


and p, = M^s\x). 

Remark 1.1. (HI) is physically motivated for classical hydrodynamics (c.f. 
mm)- {h 2) is the sound speed condition. {H3) and {H4) are the nonlin¬ 
earity conditions. 


Now, we introduce the derivatives combination of {t,u): 


y-.= y/c{u + h)^ +^ - I, q-.= y/c{u-h)^-^ + 1, (1.12) 

y c y c 


Vu 


where 


h = j ^J-pd^,Six))d^, I = 


and ho is a constant. 

Under the above assumptions, we can present the main theorem of this 
paper. 


Theorem 1.2. Suppose that (H1)-(H4) are satisfied, assume that the ini¬ 
tial entropy S = S{x) is , finite piecewise monotonic and has hounded 
total variation, if {to{x),uo{x)) are functions, and there are positive 
constants koi and ko 2 such that 


II ( to ( x ), no(x)) llci < koi, To{x) > ko 2 , for x £ M. 

Then, the solution of the Cauchy problem of ( |1.4[ ) blows up in finite time, 

if 


inf{y(x,0),g(x,0)} < -N. 

xGM. 


(1.14) 


Here N is a positive constant depending only on koi, fco 2 ond initial entropy 
function. 


This paper is organized as follows: in section 2, we introduce some nota¬ 
tions and prove the properties of the pressure p. In section 3, we obtain the 
boundedness of the Riemann invariants, and this gives the upper bound 
of density and wave speed. In section 4, we prove the finite time singularity 
formation by analysising the Riccati type equations. 
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2 Notations and Preparations 

We denote the forward and backward characteristics by 

dx dx 

— = c and — = —c, 
dt dt 

and the corresponding directional derivatives along the characteristics are 
d+ := dt + cdx and d- := dt - cdx- 
Then we can denote the Riemann invariants by 

s := u + h and r := u — h. (2-1) 

We can easily get the following system of u and h {c.f. my- 


\ht + cux = 0, 

[ut + chx +p^ = 0. 

Thus, direct calculation shows that 

d+s = {ut + chx) + {ht + cUx) = 

dS 

d-r = {ut + chx) - {ht + cUx) = 


( 2 . 2 ) 


(2.3) 


Furthermore, (1.9) yields the following inequality: 
p m'{x) 


< d+s = d-v < — 


p m! {x) 


(2.4) 


ki m{x) ^ k 2 m{x) 

Next, we will prove the property of the pressure p which will play a vital 


role in the proof of Theorem 1.2 


Lemma 2.1. Under the assumptions of (HI) and (H3), 

1 / N 1 / 

—c{s -r) <p< -c{s - r), 

where k> 1. 

Proof. If we can prove that 


(2.5) 


p < ch < kp, 
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then, due to s — r = 2h, we deduce (2.5). 

For the first part, we know that —pr is monotone decreasing in view of 


(1.5) 2 , so we have 

^+oo 

ch = 


r*+oo 




where we have used (1.11), (1.13)i and (1.5)4, thus we get p < ch. 
For the second part, if there is a constant k > 1 such that 


( 2 . 6 ) 

integrating both sides from r to +00 with respect to the time variable yields 
ch < kp. Actually, direct calculation shows that 


P'' 

C/ T 


= -\/^+ JpPrr(- 


this yields that (2.6) is equivalent to (1.8)i. Thus, 
this lemma. 


-Pt) 2 , 

we prove the result of 

□ 


Remark 2.2. It is worth noticing that p < ch is the direct conclusion of 
the property of p. So, the following property is natural 


P 

c 


< h. 


(2.7) 


Here p, c and h are defined by (1.3), (1.11) and (1.13) respectively. 


3 The boundedness of s and r 


In this section, we will first prove the L°° boundedness of the Riemann 
invariants s and r. Based on this, we can get the boundedness of |tt| and 
\h\. Finally, the upper bound of density p will be obtained, which is crucial 
to gain the singularity formation. Notice that 

d+m = cm'{x) and d-m = —cm'{x). (3-1) 


We will discuss from two aspects according to the sign of the derivative of 
m{x): 

(i) When m'{x) > 0, using (2.4) and (2.5), we have 


c m'(x) 

- (s 

2ki m(x) 


r) < 9+s 


d-r < 


c m'{x) 

2 kk 2 m{x) ^ ^ ’ 
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which means. 


1 d+m N ^ 1 9+m 

[s — r) <a+s < — 777 ^- (s — r), 


2kk2 m 


2ki m 
Introducing new variables 


2ki m 

1 d-m, s „ 1 d-TTi 

(s-r) <d-r < — — —[s-r) 


2kk2 m 


511=771^*15, rii= 777^*1 r, S12 = 777^**2 ri2 = 777^**2 r, 


then 


1 (9+777 1 (9+777 

(9+511 > 7:7 -m, (9+S12 < 7 T 7 -,-n 2 , 


2ki 777 


2kk2 777 


1 (9-777 1 9-777 

9 -rii > —-511, d-ri2 < 777^ -■S 12 . 

2kl 777 zfefe 2 ((T- 

(ii) When m![x) < 0, we similarly have 

1 (9+777 1 (9+777 

(9+S21 > 7T71-^"21, <9+522 < 7T1-^22, 


2kkl 777 


2/c2 777 


1 9-777 1 9-777 

9-r2l > 7— -521, 9-r22 < 7+; -S22, 


2kkl 777 


2/02 ((T- 


where 


521=7772**15, r2i = 7772**1 r, S22 = (^^^'“2 5, 722 =777 2*2 r. 

According to the assumptions on the initial entropy in Theorem 
have 


1 ^+00 

V:=— \S\x)\dx = 

J— 00 J—i 


7+00 

777' 

/ 

— (x) 

/ —00 

777 


dx < +00. 


Due to the assumptions on the initial data in Theorem 1.2 from m 
there exist positive constants kmi and kmr-, such that 

0 < kmi < m{x) < kmr • 

Also, we denote positive constants kg and kr, then 

|5(-,0)| </os, |r(-,0)| </Cr- 


(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 
1.2 we 

(3.8) 

: 777 ( 5 ), 

(3.9) 

(3.10) 
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Lemma 3.1. Under the assumptions of Theorem 1.2. given a point {xi,ti), 
suppose the solution of (1.4) is in the characteristic triangle bounded by 
the forward and backward characteristics through {xi,ti) o,nd the line t = 0. 
Then, we can prove that 


|r(xi,ii)| < Ur, i = 1,2. 
initial data and the number of piecewise 

backward characteristics through a point 

< X*} = {{'^{t),t)\t < t*}, 

> X*} = {(lF(t),t)|t < t*}. 

First, we will prove this lemma using the following case of three piecewise 
monotonic regions: suppose there is a point (x 2 ,t 2 ) hr the forward charac¬ 
teristic line and (xsjts) in the backward characteristic line. Assume m' < 0 
in the domain where the region from x = X 2 to x = X 3 intersects with the 
characteristic triangle, and > 0 in the rest of the characteristic triangle. 



|s(xi,tl)| < Us, 

Here Ug and Ur depend on the 
monotonic regions. 

Proof. Denote the forward and 
(x*,t*) by 

■= {(x,"?(x))|x 
■= {(x, V(x))|x 


Figure 1: Characteristic triangle 

In the forward characteristic line due to m'{x) > 0 from ("^ 1 ( 0 ), 0) 
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to (x 2 ,t 2 ), SO integrating (3.3) along this part, we can get 


1 -r 


511 ( 2 : 2 ,^ 2 ) > sii('^i(0),0) + — r —{x)rii(^x,~t{x)^dx, 

^ ^i(O) 

512 ( 2 : 2 ,^ 2 ) < si 2 C^i( 0 ), 0 ) + ^ ^{x)ri 2 (x,^{x)^ dx. 

^i(O) 

Due to the monotone increasing property of m{x), we have 

<1 for X G ( 1 ^ 1 ( 0 ), X 2 ). 


m('T’i(0)) 

< 1 

m(x) 

m{x2) 


m(x2) 


Then, we have 


m 


2: 


m 


r^x,”?(x)^ dx, (3.11) 


|5(x2,t2)| < |5C^i(0),0)| + fcs 

^i(O) 

where fca = max |( 2 A;i)“^, (2fcA;2)“^ V 

Since m'{x) < 0 , integrating ( 3 . 5 ) from (x2,t2) to (xi,ti) along the 
forward characteristic line 'ixi , we have 

52i(a:i,ti) > 52i(x2,t2) + ^^{x)r2i (^x,’?(x)^ dx, 

Xl ^ 

522(2:1, ti) < 522(2:2,^2) + ^ ^^ix)r22 (^x,”?(x)^ dx. 

X2 


Recalling (3.9), one can yield 


[m(^i(0))] 

^ kmr 

[m(x)] 

[m(xi)] 

kml 

[m(xi)] 


< 


kml 


Considering these two facts, substitute (3.7) into the above two inequalities, 
we obtain 


Xl 


|5(xi,tl)| < h\s{x2,t2)\ + h 


m' 


— (x) 

r X, t (x) 

m 

V / 


dx, 


(3.12) 


X2 
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= max 


[{2kki)-\kr, 


where k^ = max | {kmrk^i) , {kmrk^i) ^'“2 | and k^ 

From the above analyses, replacing the integration variable x by in 
(3.11) and ( 3.12[ ) shows: 

|'S(a^i,ti)| <A:4 |sC^i(0),0)| 


X\ 

+ {k^kii^ + ^ 5 ) ^ 


^i(O) 


m!, , 


- [Xa) 

r ( Xc., t [Xa) ) 

m 


dXn 


(3.13) 


As showed in Figure there are four different cases about the position of 
£ , then we will have four different results about the relationship between 

r{xa,ta) and r('fe’o-(0), 0) by the above mentioned method. When all these 
four relsuts are taken into account together, we can get 


\rixa,ta)\ < max{l,fc6}k(1r c.(0),0)| 

Xa 

+ maxl^g/cy + /cs, kj + ks} ^ 


m 


m 


(x) 


•(— 


s(x, t a{x))\dx, 


(3.14) 


aTo-lO) 


where fee = maxj {kmrk^])^'‘^ , ^7 = max{( 2 fcA;i) \ ( 2 ^ 2 ) 


and kg = max • 


:{{2ki)-^ {kmrk-])^^^^ ,{2kk2)-^ 

Finally, denoting kg = k^k^+k^, kig = max{l, fcg} and kn = max{A: 6 A; 7 + 
k%, kj + ks}, by substituting (3.14) into (3.13|), we can get 


XI 

|s(a;i,ti)| <fc4|sC^i(0),0)| + fcg/cio ^ 


Xi 


^i(O) 

Xr, 


m 


m 


r(tr^(0),0)|dxo 


+ kgkii 


m' 


m' 

- [Xa) 


— X 

m 

J 

m 


|s(x, t „{x))\dxdxo 


(3.15) 


it 1(0) 1jrCT(o) 

The first two terms can be bounded by our initial bounds. Similarly, we also 
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have 




|s(x^,t^)| <A;i2|sC^^(0),0)| + kiohs 


^e(o) 


m 


— i^c) 

m 


r(lF^(0), 0)|dX(;' 


Xf. 


(3.16) 


+ kiiki'i 


^e(o) 


m', . 

J 

m', , 

—in) 


— X 

m 

J 

lFf(0) 

m 


|s(x, t i^{x))\dxdx(^, 


where ki 2 = max{l, ^4} and ki^ = max{fc3/c4 + k^^k^ + k^}. Multiplying 
(3.16) by ^(x^), and integrating the product from xi to 1 Fi( 0) along Cxi, 
we have 


^i{0) 


m 

—i^a 

m 


Xi 


^i(O) 


< ki2 


s(x^, t^)\dx^ 

s('^g(0), 0)|(ix^ 


m 

— 

m 


Xl 


^i(O) 


+ ^10^13 


Xl 

^i(O) 


+ fcllA:i3 


m 

m 


m 

m 


X^ 


(3.17) 


^e(o) 


m 

— i^C) 

m 


|r('t^(0), 0)\dx(^dx^ 
1^c(o) 


m 

— 

m 


Xl 


Set 


^e(o) 

1 e ,,(0) 


m 


— (x) 
m 


•(— 


s(x, t (^{x))dxdx(^dx^. 


F(x,):= ^ 


m 

— X 
m 


|s(x, t rj{x))\dx. 


Since ^^(0) = ^a{0) and x^ > x^ in the same characteristic line, so we 
haveF(x() < F{xa). Combining with (3.8) and (3.10), we can rewrite (3.17) 
as 


Xl 


F{xi) < ki 2 ksV + kiokiskrV"^ + kuki^V 


1^1 (0) 


m 

— {Xa) 
m 


F{Xa)dXa. 
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Now, using the Gronwall inequality, we can get 
F{xi) < {kuksV + kioki3krV^) 

For (xo-,ta) G ~^xi, so F{xa) is also bounded by the same quantity. Thus, 
(3.15) yields that 

|s(xi,ti)| < k^ks + k'jkiokrV + k^knV {kuksV + kiokiskrV^) . 

Similarly, we can get 

\r{xi,ti)\ < kekr + kgkiok.V + kgkuV {ki 2 krV + kiokisksV^) 

From the above analyses, we can show that the Riemann invariants s 
and r are bounded in finite piecewise monotic regions. □ 


Corollary 3.2. Under the assumptions of Theorem \l.S^ we can get the 
bounds of u and h. Also, we have the upper bound of p, c andp. 


Proof. First, (2.1) gives the bounds of u and h. Due to 

-1 


-1 


ciOdf = / S)dC <h = -{s-r) 


and the assumption (1.6), there exists positive constants Tmin and Cmax de¬ 
pending only on the initial data such that 

t{x, t) > Tmin > 0, c{x, t) < Cmax- 

So we have the upper bound of p and p on account of the definition of r and 

(p^. □ 


4 Singularity formation 

First, we recall the characteristic decompositions. By the definition of y and 
q, we have {c.f. m) 

d+y = ao + aiy - a 2 y'^, d-q = ag - aiq - a 2 q^, (4.1) 


where 


flo — —cl^ + 


^ (Pr\ 

2 \ c )h 


P/i-c 



j2 

2^/c 


ai = -{2y/cl)h, 02 = 7 ^ > 0 . 

2vc 
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4.1 Estimate on the root of oq + aiy — 02 ?/^ = 0 
The first major step is to prove the lower bound on the roots of 

ao + 012/ - 022/^ = 0, 


which is 


Vroot — 


-ai ± ■\/a{ + 4ao02 _ 1 
— 2^2 2 


02 V V 02 


02 


Here (^)^ + 4^>0. It is noticeable that the lower bound of yroot depends 
on the L°° estimates of s and r. 


Lemma 4.1. Under the assumptions of Theorem 1.2, there exists a positive 


constant N only depending on the initial data such that 


1 

oi [7 

— 


2 

02 V V 


02 




02 


< N. 


Proof. We only need to show the boundedness of 0102 ^ and 0902 First, 
we give elaborative calculation on Because 


then 


0-7-T/^ C C7- Pt^ Ptt : 


{PuC ^)h = ^ ^Pfih-C ‘^ChP^i = i-Pr) - i(-pr) VrP/.- 


Also we have 


this yields ph = c, 


= -Pt = -Phhr = cph, 


Pfih = Phil = Cii = -^{-Pr) ^^Ptii, 


which implies, 


iPu^ )h~ 2 ^ 2 ^ PttPii- 


(4.2) 
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(4.3) 


Then, from the definition ( [1.13 )2 of I, we can get 

1 
4 


4 = IVc {p^^C = *PTf, - ]{-Pr) '^PttP^,- 


2 ^ /ft 4 

Integration by parts yields that 


I=- 

2 


rh 1 / 1 \ ^ 1 3 

/ Vcd{p^c~^) = - (c~^Pfj] c~2ChP^,dh 

Jho ^ ^ ^ ho ^ Jho 


{-Pt) ^Ppi 


(4.4) 


fto 


'ho 


( 

{-pT)~^PTTPfj.dh. 


Then, we have 


4 =v 


^(-Pt) + (-Pm) *P,,^, 


ho 


1 


rh 


° Jho L 

Direct calculation shows that 


7 7 T 

^(—Pt) Pt/iPttPii + {~Pt) ^PrrfiPfi T {~Pt) ^PrTPfifi 


(4.5) 


ai 


— = -2/ - 4cc^ 4 = - 2/ - S{-pr)^Prrh 

0,2 

= -2I - 2{-pr)^p7^Pri, + 2{-pr)~^Pf,. 


dh. 


(4.6) 


and 


^ = -2c^4^4 + c 4^ (PftC i)^p^-2cic^i (p^c 
= -2cic^i [4 + (p^c-i)^/] (p^c-i)^p^-4 

= - 4(-Pr)SPrr |4 “ \ [(-Pt)"Vm + (-Pr)"VrPM] ^ 
- (-Pr)^P4^PrMPM “ (-Pt)“^P^ “ 4 


(4.7) 


([ni)> 0 

^■P ^ ^ — 4^4 


By using of 

boundedness of aia^^ and ooa^Therefore, we prove this lemma. 
On the basis of the above lemma, it is easy to get: 


.9|), (|1.10|), (|2.7|) and corollary 3.2[ we can get the 

□ 


Lemma 4.2. Under the assumptions of Theorem 1.2, we ean prove that 

y{x,t)<Y, q{x,t)<Q, (4.8) 

where 


Y = max < N, sup y(x, 0) 


Q = max < N, sup q{x, 0) 


(4.9) 
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4.2 Time-dependent lower bound on 02 

To show the formation of singularity, the key step is to obtain the lower 
bound of 02 - In fact, the function 02 might vanish as time tends to infinity, 
such as for the gas dynamic case (c./. laisiiizj). 


Lemma 4.3. Assume that the pressure satisfies the assumptions (HI) and 
(H3), then 

02 {T{x,t),t) dt = 00 . (4-10) 



Proof. We know that (|4.10) is true if we can prove 


[a2(r(x,t))] ^<kut + ki5. 


(4.11) 


Direct calculation shows that 

C-h C-rTji 


02 (r(x,t)) = ^ ^ ^ (-p^) ^^Prr- 

2i/c c2^Jc 4 


Then 


[02 (r(x, t))] ^ = 4 {-Pr) 4 (prr) 


-1 


From ( 1 . 8 ) 3 , we can get 


^{-Pt)^{Ptt) <A{-Pr)*. 


(4.12) 

(4.13) 


And we also have 

( PT{x,t) 




rr(x,t) 




. t 


= \fcTt = VcUx = —{Sx+ rx) = -{y + q) 

<\{y + Q)- 


Integrating the last inequality with respect to t yields 


rT{x,t) 


rr{x,0) 


[-pdO]‘^d^< I [-pdO]‘^ d^+^{Y+Q)t < kiQt+kiY. (4.14) 


Combining (4.12), (4.13) and (4.14), we can complete the proof. 


□ 
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4.3 Singularity formation 


In this subsection, we will prove the main theorem 

Proof of Theorem 


»rem |1.2| 

is siinilar. 


We just consider the inf w(x, 0) < —N case, 

xSF 

the other case for q 

We can assume that —N is a uniform lower bound for the roots of oq + 
aiu — (1 — e)a 2 y^ = 0, according to lemma 


4.1 


Then since 02 > 0, we have 
oo + CLiy “ (1 “ < 0, for every y < —N. (4-15) 


According to the definition of infimum, there exist 0 < e <C 1 and xq G M 
such that 

y(xo,0) <-(1+ e)W (4.16) 

Now we consider the forward characteristic passing (xo,0), we have 

d+y < -ea 2 y^, 


integrating the last inequality from 0 to t with respect to the time variable, 
we can get 

t 


> 


y(xo,o) 


+ S 


^a2{l^it),t)dt, 


From (4.10) and (4.16), we can show that y blows up in finite time. 
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